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ABSTRACT 

Spectral differences between lensed quasar image components are common. Since lensing is intrin- 
sically achromatic, these differences are typically explained as the effect of either microlensing, or as 
light path time delays sampling intrinsic quasar spectral variability. Here we advance a novel third 
hypothesis: some spectral differences are due to small line-of-sight differences through quasar disk 
wind outflows. In particular, we propose that variable spectral differences seen only in component A 
of the widest separation lens SDSS J1004-I-4112 are due to differential absorption along the sightlincs. 
The absorber properties required by this hypothesis are akin to known broad absorption line (BAL) 
outflows but must have a broader, smoother velocity profile. We interpret the observed C IV emission 
line variability as further evidence for spatial fine structure transverse to the line of sight. Since out- 
flows are likely to be rotating, such absorber fine structure can consistently explain some of the UV 
and X-ray variability seen in AGN. The implications are many: (1) Spectroscopic differences in other 
lensed objects may be due to this "lens-aided multi-angle spectroscopy" (LAMAS). (2) Outflows have 
fine structure on size scales of arcsec as seen from the nucleus. (3) Assuming either broad absorption 
line region sizes proposed in recent wind models, or typically assumed continuum emission region sizes, 
LAMAS and/or variability provide broadly consistent absorber sizescale estimates of ~ 10^^ cm. (4) 
Very broad smooth absorption may be ubiquitous in quasar spectra, even when no obvious troughs 
are seen. 

Subject headings: quasars: absorption - gravitational lensing - quasars: individual (SDSS J1004-I-4112) 
- quasars: absorption lines 



1. INTRODUCTION 
1.1. Quasar Outflows 

Evidence is accumulating that outflows occur wher- 
ever there is accretion. In active galactic nuclei (AGN) 
the outflows from supermassive black holes (SMBH) are 
highly ionized, so their signatures appear mostly in rest- 
frame ultraviolet (UV) and X-ray spectra. From spec- 
troscopy of Seyfert nuclei, at least half show narrow ab- 
sorption lines (NALs) of highly ionized species in out- 
flows of ~1000 km s^^ . From X-ray studies, warm ab- 
sorber feature s are present again in about half of quasar 
X-ray spectra (jPiconcelli et al.l2005D. and many such fea- 
tures are seen in outflow jCrenshawe t al. 2003a). 

More massive, accelerating outflows create broad ab- 
sorption lines (BALs), whose P-Cygni profiles span 
velocities to ~0.3c and are visible in the spectra of 
15 - 20% of o ptically-sele cted quasars ()Hewett fc Foltd 
l2003t iR eichar d et al.li200^ or more in unbiased samples 
l)Chartasll2000l|) . These outflows may exist in all quasars, 
subtending a solid angle covering fraction at least as large 
as their detection fraction. 

Studies of powerful mass outflows in quasars are 
rapidly reshaping our understanding of physics near the 
supermassive black hole. Sightlines that pierce the ab- 
sorbers yield information on the velocity and plasma 
state of the gas that is impossible to obtain from emis- 
sion lines, most likely because the latter are formed by 
compounded emission from much large r ensembles of 
clouds in a wide variety of physical states l|Baldwin et al.l 
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Il995|l . These broad line-emitting clouds lie at distances 
of 10^^ — 10^^ cm from the SMBH, spanning a range of 
densities ~ 10* — 10^^ cm~^ and covering >10% of 
the ionizing source. Recent models (e.g., ;;Elvis 2000) 
hold that the absorbing and emitting clouds are not only 
cospatial but quite possibly identical, making knowledge 
of absorber properties even more important to our un- 
derstanding of AGN physics. 

For outflows, equato rial disk wind models (e.g., 
iMurrav fc Chiand Il997t IeI^ |2000|) are currently fa- 
vored, although proponents of polar wi nds exist 
llPunslv Lina,rill2nnf iH^uoll BlackmanI 120011). It 
has often been suggested that outflows may con- 
tain a substantia l kinetic luminosity (though see 
iBlustin et al.ll2b0!Tfl , and impart significant mass and en- 
ergy into the interstellar (ISM) of the host galaxy and 
even into the surrounding in tergalactic medium (IGM; 
iRovchowdhurv fc Nathll200^ . 

A great deal more information about the structure and 
dynamics of the outflowing absorber could be gleaned if 
some information about their size scales were available. 
However, the ~parsec scale emitting/absorbing regions 
near supermassive black holes (SMBHs) will remain spa- 
tially unresolved for the foreseeable future (milliarcsec at 
Z'^O.Ol). If the dynamics of absorbers were known (e.g., 
Keplerian orbits in the SMBH potential), then spectro- 
scopic variability information might also lead to absorber 
size scales, viz. ay ~ Vtrans^t- Unfortunately, the ab- 
sorber dynamics are poorly constrained, and multi-epoch 
spectroscopy is difficult to arrange and therefore rare. 
Gravitational lensing can help. 
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1.2. Spectral Differences Between Lensed Images 

Because gravitational lensing is intrinsically achro- 
matic, spectral similarity is an early criterion to con- 
sider close quasar images as lens candidates. Redsliifts, 
as well as emission and absorption line profiles must be 
similar. However, significant spectral differences have 
been noted in bona fide l ensed quasars (e. g., HE 2149- 
2745; Burud ct al. 2002a, SBS 1520-1-530: Buru d et all 
l2002b: Ogu ri et al. 2005) with clearly identified lensing 
masses and some with time delays (making the lens in- 
terpretation irrefutable). Significant differences in ab- 
sorber properties between image components have been 
documented in optical/ UV spectra of BALQSO lenses 
(e.g., APM082 9-f5255: ILewls et alJ .2002. H1413-H17; 
lAngonin et al][l 990'). SDSS J1004-h4112 could represent 
a more extreme example of this phenomenon due to its 
wider- angle separation, and/or to fortuitous snapshots of 
unveiled phases. 

One traditional explanation of such spectral differences 
is microlensing, which preferentially magnifies parts of 
one image, enhancing spectral components that originate 
from smaller emitting regions (about the size of the Ein- 
stein radius of a star) at the source. Another possible 
explanation is intrinsic quasar spectral variability, com- 
bined with lensing time delays. 

The recent discovery of the gravitationally lensed 
quadruple-iinage z = 1.734 quasar SDSS J1004+4112 
l|Inada et al.l l2003': T) gu ri et aDl2004D is exciting because 
of its record-setting separation ( maximum of 14.6" 
between images). A z=0.68 cluster centered among 
the four lensed images is confirmed as the massive 
lens. More recently, deep ACS and NICMOS images of 
SDSSJ1004+4112 from Hubble Space Telescope (HST) 
(jinada ct al. 2005) have now revealed clear arcs, sheared 
images of the quasar host galaxy, and a probable 5th 
quasar image, all of which substantially constrain the 
plethora of viable lens models. Here we discuss primar- 
ily the 4 brighter image components A — D with existing 
spectroscopy. 

In addition, there are intriguing differences between 
the spectra of the 4 quasar images. Both microlensing 
and variability have been posited to explain the spectral 
differences, yet both explanations have serious problems. 

Here we advance a novel third hypot hesis that has 
only been mentioned in passing (e.g., iLewis fc Bellel 
ll998HOguri et al.l2004fl - line-of-sight differences through 
quasar outflows. We propose that SDSS J1004-t-4112 of- 
fers a revealing multi-angle view of quasar winds origi- 
nating near the nucleus, where absorption can change on 
small angular and time scales. While there are challenges 
also for this new hypothesis, several of its predictions are 
immediately testable. Given the wide ramifications for 
AGN physics, it is well worth considering. 

In the next section, we briefly review the spectral ev- 
idence, and troubles with the more traditional explana- 
tions. Then we examine the prospects and predictions of 
our own hypothesis. 

2. SPECTRA OF SDSS J1004+4112 

2.1. Description 

In SDSS J1004-I-4112, the blue emission line wings 
of the brightest image A are enhanced relative to the 
spectra of image components B, C and D. The ra- 



tios of compone nt spectra reproduced in Figure ^ 
from iQguri et al.l (j2p04) reveal that the differences are 
larger for LyaA 1216/N VA 1240, SiIV-HOIV]A 1400 and 
CIVA1549 lines than for the lower ionization lines of 
CIII]A 1909 and MgllA 2800. There are subtle differences 
between the ratio spectra A/B, A/C, and A/D, but the 
general features are similar; the spectrum of A is the 
most divergent, with strong blue wings to the emission 
lines. ^ The continuum flux ratios are approximately 
constant from 3000 - 8OOOA. ^ 

Figure direct from from iRichards et alJ l)2004j) . 
focuses on multi-epoch spectra of the CIV emission 
line in components A and B. The 7 spectra shown 
span observed-frame time delays of 322 days. An en- 
hancement of the blue wing of C IV in A was seen 
in the first (03 May 2003) spectrum which lasted at 
least 28d (since it is seen on 31 May), and then faded 
(since it was not there 21 Nov). Not shown in the 
figure are later epochs when the enhancement reap- 
peared — seen in a spectrum of 26 Mar 2004, and again 
10 and 19 April 2004 (,Richards. Johnston, fc Hennawi 
l200llWisotzki et al.ll200ir 



2.2. The Microlensing Hypothesis 

Could the blue wing enhancement in A be due to mi- 
crolensing? Microlensing of the broad emission line re- 
gion (BELR) can occur if it has structure smaller than 
the Einstein radius (~ 3x 10^'"' cm for a O.IM0 star). Sev- 
eral problems plague the microlensing explanation: [1] 
Microlensing should amplify not just the BELR, but also 
the hot continuum-emitting region interior to it. While 
the single star approximation is poor for microlensing at 
significant optical depth (the situation for multiply im- 
aged quasars) , caustic models show a strong general cor- 
relation between the magnification of the continuum and 
of the BLR, and it is extremely r a.re for the BLR to b e 
magnified but not the continuum l)Lewis fc Ibatal l200l. 
No amplification of the A continuum was seen (<20%; 
Richards et al. 2004) and t he A and B continuua are 
effectively identical ( Wisotzki et al.ll2004(l . [2] If the mi- 
crolensing hypothesis is correct, then microlensing does 
not act on the (smaller) continuum region but somehow 
acts only on a select region of the BELR, and that this 
same configuration recurs. The reappearance of the same 
enhancement renders the microlensing explanation par- 
ticularly unlikely. [3] Strong line profile differences are 
also seen in the blue wings of the lower ionization lines 
of CHI] and Mgll (Richards ct al. 2004). In most BELR 
models, these come from significantly larger regions than 
does C IV emission, and should be less susceptible to mi- 
crolensing. [4] Line asymmetry induced by microlensing 
is expected for certain kinematic models of the BELR 
(Abaias et al. 2002), the closest match being for a mod- 
ified Keplerian BELR with small lens size. The A profile 
has varied strongly. Similar profile variations would be 
expected in other lensed systems but have never been 

^ We not e that , while not relevant to the present discussion, 
IQguri et alj 120041) discuss also the variable strengths of narrow 
intervening e.g., Mgll absorbers that can be seen in these spectra. 

^ Not all spectra were taken at the parallactic angle, so the 
continuum slopes are suspect. However, neither would we expect 
these parallactic angle effects to conspire to produce flat continuum 
ratios. 
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seen.^ 

We note that in the ACS and NICMOS images 
ijinada et al.ll2005j) . a dim object (G4) is detected within 
^1" of image A that could host a microlens. However, 
it is uncertain whether or not the object is a galaxy, and 
even more so whether it could provide a microlensing op- 
tical depth at the position of image A sufficient to explain 
the observed variations. 

2.3. Variability 

Spectral differences in many lenses are plausibly ex- 
plained by intrinsic variability combined with time de- 
lays, meaning that with 4 images, we are effectively 
viewing one quasar at 4 epochs. Given the asymmetry 
{rA — tb)/ [ta +tb) of the A and B images with respect 
to the lens, the maximum delay between them is <, 30 
days l|Oguri et al.l f2004). However, B never showed a 
blue wing bump, althovigh it persisted in A for at least 
a month (Richards ct al. 2004). Also, the persistently 
bluer profiles in A should have appeared within ^^BOd in 
B (or disappeared in A depending on which image lags) 
but have not: the variability has been confined to A. For 
these reasons, the usual hypothesis of intrinsic variability 
plus lightpath time delay appears unlikely to be correct. 
However, the temporal coverage of spectroscopy to date 
can not e ntirely rule out intri nsic variability plus time 
delay fsee iRichards et aDl2004l for a full discussion). 

2.4. Lens-Aided Multi-Angle Spectroscopy through 
Absorbing Fine Structure 

Lensing geometry actually provides slightly different 
sightlines to the lensed object, and the image sepa- 
ration 9 is similar to the angular difference between 
the lines of sight as seen from t he quasar nucleus 
([Schneider. Ehlers. fc Falcd IT992(l . Indeed. IChelouchd 
()2003|) recentlv pointed out the utility of lensed images 
of BALQSOs for probing absorber structure transverse 
to the line of sight. 

We propose that small angular differences in sight- 
line afforded by lensing (lens-aided multi-angle spec- 
troscopy; LAMAS hereafter) can probe significantly dif- 
ferent absorbing columns in quasars. In the case of 
SDSS J1004+4112, we propose that all 4 image com- 
ponents suffer from absorption from a warm outflowing 
wind (similar to that proposed by Elvis 2000 and others), 
and that the spectral differences are due to a line-of-sight 
to component A that pierces a persistently thinner and 
patchier part of the absorbing flow. A transient hole or 
"rend" in the flow along that sightline is responsible for 
the blue wing C IV enhancement in A. 

We note that a surprising Galactic analogy to LAMAS 
exists in the massive young star r) Carinae. There is 
a dusty axisymmetric bipolar "Homunculus" , formed 
by ejecta from r] Car, which creates a hollow reflection 
nebula. Spatially -resolved HST STIS spectra of the 
Homunculus (Smith et al.l l2003j) , whose 3-dimensional 
structure is fairly well-modeled, samples reflected light 
from the star emerging at different stellar latitudes and 
probing different parts of the outflow. At viewing angles 
separated on arc sec scales, large differences in the broad 
(~I000 km s~^) P Cygni absorption proflles are seen. 

3 While SDSS J1004-I-4112 is uniquely wide, this is irrelevant to 
microlensing. 



Outflow structure in a young massive star may well be 
qualitatively different than in an AGN, but the dense 
line-driven wind and the strong angular dependence of 
the absorption are analogous. 

First we present evidence that the blue enhancements 
in SDSS J1004-h4112 A could be due to the alleviation 
of absorption along this sightline. Second, we report on 
significant variations common in AGN absorbers. Third, 
we address the geometry - how different lines of sight 
may probe significantly different absorber properties and 
what constraints this and variability place on the ab- 
sorber structure. Fourth, we address some possible ob- 
jections to our hypothesis. 

3. RENDING OF THE VEIL 

The ratio spectra shown in Figure^ are strongly rem- 
iniscent of the structure seen in BALQSOs. To illustrate 
this with typical spectra, we take the the rest-fram e non- 
BAL composite spectrum from the SDSS ( Rcichard et alJ 
120031) . assumed to be "unabsorbed" . Next we gener- 
ate "part-BAL" spectra by combining that composite 
with the SDSS high-ionization BAL composite spectrum. 
The non-B AL spectrum is shown at the top of Figure 13 
along with three example part-BAL composite spectra 
generated by adding 25, 20, and 15% of the SDSS high- 
ionization BAL composite spectrum. Even though we 
plot the logarithm to make differences as visible as pos- 
sible, no clear absorption signature is evident even in 
the 25% part-BAL spectrum. Then we divided the part- 
BAL composite by the "unabsorbed" (non-BAL) com- 
posite (lower panel of Figure [2Jl . While some features 
may differ in detail, note the overall similarity to the 
spectral ratios seen in Figure[l]for SDSS J1004-I-4112. 

The BAL and non-BAL composites have all their fea- 
tures smoothed and broadened by the averaging process. 
These composite ratio spectra (of a part-BAL to a non- 
BAL composite) in Figure |21 show strong similarities to 
the SDSS JI004-I-4112 image component spectral ratios 
in Figure ^ This implies that absorption may account 
for the differences along the sightlines to the different im- 
ages. Furthermore, there are 2 reasons to think that the 
absorbing outflows in SDSSJ1004 may be smoother and 
broader than normally associated with BALs: (1) adding 
a smooth/broad absorbed spectral component as in Fig- 
ure|21(T0P) does not create recognizable absorption sig- 
natures, nor are such signatures visible in the BCD spec- 
tra. (2) The similarity of the ratio spectra in Figures 1 
and 2 suggest that the absorption troughs in the BCD 
components may also be smoothed and broadened, as 
are those in the composites. Normal BAL-type absorbers 
typically show steep trough edges, which are clearly ab- 
sent in SDSS J1004+4112. While the absorption in the 
BAL composite is smoothed by averaging many sightlines 
to different normal BALQSOs, the absorber we hypoth- 
esize to exist in SDSS J1004-I-4112 must be smoothed 
in velocity for a different reason. We propose that the 
absorber profile is intrinsically smeared broadly in ve- 
locity and consequently smooth onset and trailoff, i.e., 
with a profile more like a shallow bowl than a trough. 
Because such absorption is easily missed, this may be a 
higher covering-factor but lower-column counterpart to 
more typically remarked BAL flows. 

We propose that virtually all quasar sightlines pene- 
trate a warm (ionized), relativistically outflowing wind 
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which is smooth in velocity space but may be spa- 
tially clumped or filamentary. These ionized winds 
would quite normally sculpt the emission line profiles of 
all quasars, particularly in the UV. Excess emission in 
the blue wing flux of resonance lines is not uncommon 
llBachev et alJ 120041) . and is predicted in the models of 
iMurrav fc Chiand l|1997|l (see their Fig 7). The latter 
authors note that "observed profiles are even more an- 
tisymmetric than the calculated profiles. Our neglect of 
scattered line photons may be partially responsible for 
this. The absorption troughs in BALs provide direct ev- 
idence for the existence of such scattered line photons." 

In the case of SDSS J1004-H4112, our line of sight to 
A, as influenced by the lensing geometry, probes through 
part of the absorber whose column density is somewhat 
less dense, and also more susceptible to the appearance 
of gaps in the patchy or perhaps filamentary absorber. 

As a microlensing effect, similar blue wing enhance- 
ments would not be expected in unlensed quasars. How- 
ever as either a variability or a line-of-sight effect, they 
should be seen (with as yet unknown frequency) in the 
general (unlensed) quasar population. We investigated 
the emission line properties of AGN from the large sam- 
ple of HST FOS line measurements of Kuraszkiewicz et 
al. (2002, 2004). While their spectral analysis was both 
careful and uniform, the sample is heterogeneous, and 
our perusal subjective, so we cannot draw firm conclu- 
sions on the fraction of objects with blue asymmetry 
at this strength. However, we find several examples of 
C IV lines with asymmetry simflar the SDSS J1004-f 4112 
spectrum seen in Figure ^ The spectrum and model fits 
to PKS J2355-3357 are shown in Figure |3 This sug- 
gests that the same phenomenon occurs in other objects 
along unlensed sightlines as expected. Because we have 
no multi-epoch spectroscopy of these objects, we do not 
know whether the phenomenon is typically transitory. 

The absorbing structure we propose is not identical to 
a BAL fiow as typically observed: its velocity profile is 
broader and smoother. This is because no steep trough 
edges are seen in Figure ^ (Of course, no steep edges 
are seen in the comparison composite spectrum of Fig- 
ure 121 either, but this is due to the fact that the SDSS 
BAL composite is an average of 180 SDSS BALQSOs, 
whose BAL troughs have a range of detachment veloci- 
ties and velocity profiles.) Broad smooth outflow velocity 
profiles are plausible because BALQSOs with unusually 
wide, smooth absorption have b een found. Extreme ex - 
amples like VPMS J1342-H2840 l|Meusinger et aLll20(il 
SPSS J0105-0033 and SDSS .12204-^0031 iHall etall 
120021) lack any evidence for the usual restframe UV emis- 
sion lines, yet they show relatively blue continua with no 
obvious dust reddening. The most convincing explana- 
tion to date of their spectral features is unusually wide, 
overlapping low- ionization BAL troughs (loB ALs) . Since 
these objects are quite diflicult to recognize and clas- 
sify, they are likely severely under-represented in exist- 
ing AGN samples. From our crude spectral arithmetic in 
Figure|21and from the fact that most quasars do not show 
X-ray absorption signatures as strong as in BALQSOs, 
the flow we propose probably has lower overall column 
than typical BALs. Indeed, there is evidence that the 
frequenc y of qu asar absorption increases towards lower 
columns (jReichard e t al. 2003), making more plausible a 
ubiquitous low-column smooth absorber. 



The disk wind scenario iMurrav fc ChiandlToOTHEM^ 
120001) fits naturally into this picture, since the absorbers 
are smoothly accelerating outfiowing sheets of warm ion- 
ized plasma. In the disk wind scenario, radiative acceler- 
ation is in the radial direction, but the disk from which 
the wind arises is rotating at approxima tely Keplerian 
speed until the last marginally stable orbit l)Murrav et alJ 
Il995j) . Put simply, the outflows are almost certainly ro- 
tating, creating helical streamlines. These are thought to 
be rising off the disk to create a "martini glass" funnel- 
shaped outflow (,Elvia.200Q,') . 

We might expect that holes or fiow gaps would ap- 
pear more commonly along the sightline to component 
A of SDSS J1004+4112 if it is less absorbed generally 
(see Figure 0J. Even after the obvious blue wing bump 
subsided during later 2003 epochs shown in Figure ^p, 
the blue/red fiux ratios (measured 20 A to either side of 
CIVA1549) are Ra = 1.15 and Ra = 0.84 (both ±0.05). 
This supports our hypothesis of persistently lower col- 
umn along the sightline to the A image component. 

3.1. Variable Absorption 

Time-dependent calculations iProga et al.ll200(Tl) show 
that disk wind instabilities result in a filamentary sub- 
structure to the flow, so we naturally expect column den- 
sity variations to occur along our sightline, which would 
more strongly affect the blue line wings. Rotation com- 
bined with clumpy outflows predicts absorber variability. 
This prediction should be borne out in absorbers seen in 
unlensed quasars. 

In the UV, BALs have been seen to vary strongly 
in QSOs. At least one BALQSO spec troscopic mon - 
itoring campaign has been performed ijBarlowl Il994() . 
which showed BAL variability in 15 of 23 BALQSOs. 
Significant BAL variability has been reported anecdo- 
tally el sewher e (CSO203; Barlow et aia9_92, SDSS 0437- 
0045; HaUetal. 2002). The BALs in at le ast one 
quasar have disappeared completely l|Mal l2003) . Vari- 
ations may be due to a combination of factors, in- 
cluding a change in line-of-sight column, covering fac- 
tor of the absorber, or ionization. Variability is well- 
documented and quite common in the narrower absorb- 
ing systems seen in lower luminosi ty AGN (Seyferts) 
such as s uch as NGC351 6 ( Koratk ar et al.l [l996'l and 
NGC4151 (Kraemer ct al."'2b05) and others, as reviewed 
by Crenshaw ct al. (2003a). 

In the highly ionized circumnuclear environment, X- 
ray spectroscopy is a sensitive probe of column den- 
sity because metals in virtually any ionization state 
will absorb X-rays, whereas optical/UV continuum ab- 
sorption requires the presence of dust. For instance, 
while the continua of high-ioni zation B ALQSO s ar e 
at most only mildly reddened ijReichard et al.,, .20031). 
BAL QSOs are strongly absorbed in X-ravs ijGreen et alJ 
I200H iGallagher et al.l l2002a() . As measured by X-ray 
spectroscopy, variations in the absorbing column den- 
sity by factors of a few are quite common in both 
optical broad and narrow line (Type I and II) AGN. 
In Seyfert 2 galaxies, 20 - 80% variability in the 
measured absorbing column is endemic (23 of 24 ob- 
jects; Risaliti. Elvis. & Nicastro 2002), and is likely due 
to bulk motion of n iaterial across the line-of-sight. 
IGallagher et alJ l)2004|) found a column variation of 
6x10^2 atoms cm'^ in BALQSO PG2112-h059 over 
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a ~3yr timespan. iGallaeher et alJ l)2002bfl discovered 
hard-band variability at the 45% level on a timescale of 
20ksec in the nearby mini-BAL QSO Mkn 231. More 
dramatic changes are also seen. UGC 4203 underwent a 
transition in its X-r ay spectrum from Compton-thick to 
thin ( Guainazzi et al 2002) . and NGC3227 did the op- 
posite l|Lamer et al.ll2003iK The authors suggested the 
transition could be absorbing clouds or streams crossing 
our line-of-sight. An X-ray "unveil ing event" was noted 
in NGC4388 by ijElvis et al.ll2004|) . corresponding to a 
decrease in column of a factor of 100 in just 4 hours. 
Variability so rapid puts the NGC4388 absorber at a few 
100 i?s (~ 3 X 10^5 cm for Mbh = 10* Mq), similar to 
the broad emission line region or smaller. 

Absorber variability is clearly common if not endemic 
even in the general (unlensed) AGN population. Such 
variability should provide a geometric constraint ay on 
the absorber sizescale from A t that is complementary to 
the LAMAS sizescale from 9. 

4. GEOMETRY 

From the LAMAS perspective, the different absorber 
properties of the two sightlines constrain the lateral 
sizescale a of the absorber as simply as ^ RaQ, where 
Ra is the absorber distance and 9 is the observed image 
splitting wherein a signifcant absorbing column change 
is seen. The difSculty here is twofold. First, there are at 
best only loose constraints on location of the broad ab- 
sorption line region (BALR): debate over Rr , still span s 
5 orders of m agnitu de from 0.01 to 1000 pc <|ElviiJl200'fll: 
IdeKool et all 120011: Is^erett et all I2002M Second, the 
smooth, broad absorbing flow we propose may span a 
different spatial regime than the dense flows responsible 
for BALs (the BALR). 

We first assume that Ra ^ Rbelr- From reverber- 
ation mapping (Peterson 1997), the size of the high- 
ionization (CIV) BELR in NGC 5548 is R ~10 It-day, 
or ~ 10"'^^ cm. However, in higher luminosity objects 
like SDSS J1004-I-4112, reverberation time lags can eas- 
ily reach 100 It-day and more. (The BELR size scales 

1/2 

as Rbelr oc 0.01 L44 pc from the central continuum 
source, where Laa is the 0.1 - Ifxm luminosity in units 
of lO'^'^erg s-i (jNetzer fc PetersonllTgg?: Peterson et al.l 
|2004).) If the absorber is at a distance similar to the 
BELR (assuming 100 It-days), then the A/B spectral dif- 
ferences seen in SDSS J1004-I-4112 across 9 = 3.7" imply 
that cloud/stream columns can change significantly on 
size-scales of ~ 5 x 10^^ cm or < lA.U. transverse to the 
line-of-sight. 

If we assume that the LAMAS and variability both 
probe the same absorber size scale, how do these sizes 
compare? Suppose the flow to be in quasi-Keplerian ro- 
tation at ~ lO'' km s~^ (t he median FWHM of broa d 
component of CIV emission: iKuraszkiewicz et al.ll200'^ . 
The observed variability timescale of the C IV blue wing 
thus potentially presents a separate (A0) size constraint 
complementary to that of LAMAS {9) discussed above: 
the blue wing bum p disappearing from A in < 43days 
(in the rest frame; iRichards et al J 120041) constrains the 
transverse size of a rend or flow gap to av<,4 x lO^^cm 
(~270 AU). If we again assume transverse cloud motion 
created the variable line profile, but instead of postulat- 
ing a privileged sightline to A we make the (unlikely) 
assumption that with adequate monitoring, the same 



change could have been seen to occur in components 
B ~ D, then the lensing delay timescale is a more ap- 
propriate compariso n. Since mo st lensing models for 
SDSS J1004-f 4112 in lQguri et all (|2004) predict a time 
delay < 30 days, the result for ay is therefore similar. 
Either way, ay is inconsistent with (lOOOx) the LAMAS 
sizescale ul based on Ra ~ Rbelr above. 

Using a single-phase wind to model the outflows from 
FIRST J104459. 6+365605, dcKool et al. (2001) found 
the low-ionization BAL region must be at Ra ~700pc. 
In this case, the lateral absorber size implied by LAMAS 
is ai ~ 4 X 10^^ cm, which i s lOx larger than the variabil- 
ity sizescale. However, the IdeKool et al.l l|2001l) sizescale 
is vary hard to reconcile w ith partial covering seen in 
BALs. lEverett et al.l (|2002) instead used a multi-phase 
outflow model that successfully reproduced the many ab- 
sorption features (having different ionization parameters 
but similar velocity structure) by placing the absorber at 
^ 4pc from the central source. Using this value for Ra, 
LAMAS yields (ol 4 x lO^'^cm or 30AU), about lOx 
smaller than the variability estimate. 

So if our LAMAS proposal holds true, the absorber dis- 
tance is inconsistent with most BELR distance estimates 
of ^ 10^^ cm, but falls somewhere between the best re- 
cent estimates of the absorber distance Ra for BALs of 
5< Ra <700pc. 

As an independent estimate, we consider that from an 
absorbed sightline, the covering fraction of the absorber 
must be reasonably large (between 0.1 and 1) for a small 
change in viewing angle to significantly affect the ab- 
sorption profile. Put differently, the size of the contin- 
uum region as seen from the absorber distance must not 
be too large, otherwise a small change in viewing angle 
would not detectably alter the absorption profile. For 
SDSS J1004+4112 (with its 3.7" splitting), this implies 
that Ra > 50000 ac, where ac is the projected UV con- 
tinuum emitting region size. Because ac is thought to 
be about ~ 30 — 50 i?s, for an accreting supermassive 
black hole of 10* Mq we expect that ac ~ 10"'^^, yielding 
i?a > 5 X lO^^cm or >10pc. While independent from 
the above geometrical arguments, this covering fraction 
argument under the LAMAS hypothesis yields a simi- 
lar plausible absorber sizescale (10^^ cm). Note that if 
spectra were to show the absorptionn of both contin- 
uum and substantial BELR flux, then the absorber dis- 
tance implied by this latter argument would be quite 
large (>,1 kpc). 

5. OBJECTIONS 

Several objections to LAMAS leap to mind. 

How could it be that the absorption changes on arcsec- 
ond sightline scales but is identical in components B — D, 
some of which are more widely separated from each other 
than from A? Our proposal is that our sightline to A 
is unique in this system, perhaps skirting the edge of 
a larger-scale structure like a disk wind, as illustrated 
crudely in Figure^ However, there is no reason in the 
LAMAS picture that only one component could show a 
distinct absorption profile. Spectroscopic monitoring of 
this and other lenses may well identify such cases, which 
would help confirm the model. Our absorption interpre- 
tation would be most sensitively verified by simultane- 
ous UV and X-ray monitoring of spatially-resolved image 
components in this and other lenses. 
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SDSSJIOO4+4II2 IS not a BALQSO. Most 
or possibly all QSOs contain BAL-type outflows 
(jHamann fc Ferlandlll993j) . But the classic deep BAL 
troughs are ob served only for sightlines traversing dense 
BAL streams (Ogle' "1999), whereas lower-column parts 
of the flow cover more solid angle and may aff ect the 
emiss ioi i lines less spec t acular ly in most QSOs ijGreenl 
Il998|) . iReichard et all l|2003j) noted that the fraction 
of quasars with BALs increases strongly towards lower 
BALnicity, so that "the fraction of quasars with intrinsic 
outflows may be significantly underestimated." Broader, 
smoother outflows such as proposed here be even harder 
to detect, but could be ubiquitous. 

If such smooth, broad outflows are ubiquitous, what 
special conditions yield more typical BALs, which are de- 
tached blueward of the emission line peak velocity? Since 
their measured X-ray column densities are highest of all 
types of QSOs, the densest part of the wind is prob ed by 
BAL-t ype outflows. In a disk wind structure ala lElvi j 
(|2p00') , the arch of the BAL wind (as it moves from 
predominantly vertical to radial velocity) means that it 
is already accelerated when it first crosses our sightline 
(Figure Hi). This need not be true for the WHIM that 
transports the BAL material. The WHIM can pervade 
a much larger opening angle, and be accelerated from a 
lower velocity. 

// smooth, broad outflows are ubiquitous, why aren't 
most quasars absorbed by columns as large as '^10^^ 
in the X-ray regime? X-ray spectral fitting to ab- 
sorption features in bright optical- and radio-selected 
quasars typically yields measured columns of '^10^^ 
ijReeves fc Turneil 120001) . However, very broad fea- 
tures from warm ionized gas would not be easily de- 
tected, just as they are unrecognized in the major- 
ity of UV spectra. In fact, for Seyfert Is and radio- 
quiet quasars lacking obvious strong absorption signa- 
tures in the UV and X-ray regimes, there is a ubiq- 
uitous and poorly expl ained broad soft spectra l ex- 
cess below about 2 keV. iGierlihski fc Done! l)2004|) and 
iSobolewska fc Don^ 1)20051) have interpreted the soft ex- 
cess as an artifact of previously unrecognized, relativis- 
tically smeared, partially ionized absorption; strong very 
broad O VII, O VIII, and Fe absorption features at 0.7 - 
0.9 keV can lead to an apparent upward curvature below 
these energies, mimicking soft excess emission. Whether 
such fiows are related is still unclear. Their proposed X- 
ray absorber velocities [v/c ^ 0.2) greatly exceed what 
we see in the UV spectra of SDSSJ10 04 (vie ^ 0.04) 
but th e absorbing zo nes mav be disjoint ([Crenshaw et al.l 
l2008bD or stratified l|Steenbrugge et al.ll200,'TD . or the ve- 
locity range may vary between objects. 

Some significant details of the part-BAL com- 
posite spectra in Figure differ from those of 
SDSSJIOO4-I-4II2. Figure El is meant to be illustrative 
only, because there are several reasons why those spec- 
tra differ in some of the details: 1) The BAL compos- 
ites are rather heavily smoothed (by dint of averaging 
many BALs with different velocity profiles). Weaker fea- 
tures are thereby blended away, and then further diluted 
by the weighting in the part-BAL composite. 2) De- 
tails of the composites depend on their method of con- 
struction: how the mean is weighted, in what rest wave- 
length region the contributing spectra are normalized, 
the dereddening algorithm, and how the quasars' sys- 



temic velocity is defined for deredshifting. Furthermore, 
the SDSS sample selection and spectroscopy means that 
QSOs contributing to the composite have varying lumi- 
nosity, redshift, and intrinsi c reddening as a function of 
rest wavelength (Willott 2005). 3) The SDSS composites 
match SDSSJ1004 in neither luminosity nor redshift, so 
that related emissi on line differences are expec ted (e.g., 
the Baldwin effect lYip et al]l200l IGreenllTool in CHI] 
and Hell. 4) True BALQSOs typically show other differ- 
ences in emission line profiles and strengths relative to 
non-BALQSOs that are possibly rel ated to an enhanced 
accretion rate (e.g., iBorosonI I2001 . 5) The LAMAS 
outfiow may sample different regions than typical BAL 
fiows, perhaps with different launch point and/or accel- 
eration profile. Assuming a (e.g., biconical) disk- wind 
outfiow model, our sightline to SDSS J1004-|-4112v4 may 
be somewhat farther from the denser part of the outfiow, 
especially given the claim of repeated 'rending' of the 
fiow. In this part of the outfiow (e.g., the lower column 
warm highly ionized medium - the 'WHIM') the pres- 
sure/ionization balance may be different than in a typical 
BAL fiow. Given so many caveats, the accordance be- 
tween the BAL composites and FigureQseems rather re- 
markable. And certainly no plausible, testable combina- 
tion of emission line production and microlensing theory 
is yet available that can reproduce the SDSS J1004-I-4112 
spectra as well. 



6. CAVEATS 

Whether due to different effective orientations (caused 
by lensing-induced sightline difference) or to time vari- 
ability, differential absorption by itself may not be able 
to produce three key features seen in the spectral differ- 
ences between components. First is the apparent absorp- 
tion in CIII] A 1909 and He HA 1640. Second, the LAMAS 
model implies that emission lines in most quasars would 
have an intrinsic blueward asymmetry except for typi- 
cal broad smooth absorption caused by warm outfiowing 
winds. If true, emission-line profiles not subject to res- 
onant absorption (e.g., CIIIJA1909) should retain that 
asymmetry in the general AGN population. Finally, we 
discuss the apparent observed change in the red wings of 
the emission-line profiles (between images and also tem- 
porally). 

1) The blue wing of the C HI] emission profile also ap- 
pears to be enhanced in image A relative t o the other 
images (see Figure 4 of Richar ds et al.ll200^ , but there 
should be no absorption (and so no differential absorp- 
tion) detectable from the semi-forbidden CIII] transi- 
tion. However, this spectral difference may still be due 
to absorption from another species, since BALQSOs can 
show some spectral signatures there. This is evident 
for instance in the plots like Fig 10 of iRichards et all 
(2002), where their BALQSO composites show some sig- 
nificant C HI] line strength and profile differences from 
the non-BAL comp osite. The tro ugh seen near restframe 
1909A in Richard s et alJ l|2004D Fig 4 is not necessar- 
ily from CIII]. Felll absorption (UV 34,UV48 ) are strong 
iron lines in the CIII] region (Graha m et alJ ll996). Such 
lines hav e been see n in a bsorption in BALQSOs e.g.. 
Fig 22 of lHall et all 1)20021) . though those tend to be in 
the more rare low-ionization BALQSOs (loBAL QSOs) 
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only.'* For unknown reasons, strong Fell and Felll emis- 
sion blends are endemic to BALQSOs generally, an d cor- 
relate with BAL strength iWevmann et alJlT99H) .^ So, 
some differences in iron blend profiles are not unexpected 
if the BCD component spectra are more BALQSO-like. 

Also troubling is that the SDSS J1004+4112 ra- 
tio spectra appear to show an absorption feature in 
HellA 1640 which is rare in BALs. Some effect related to 
absorption may be at work here. For instance, for rea- 
sons as yet unknown, the He II emission feature in BALs, 
and in objects with larger CIV blueshifts is weaker and 
broader in general (c.f. BALQSOs and composite C in 
Figures 4 and 10 of lRichards et alJl2002D . 

2) In the general AGN population, emission-line pro- 
files that are not subject to resonant absorption (e.g., 
CIIIJA1909) should retain the blueward asymmetry im- 
plied by the differential absorption hypothesis. While 
this is a cogent objection to the differential absorption 
hypothesis, it has only weak empirical evidence avail- 
able to test it because there are very few strong non- 
resonant emission lines seen in AGN spectra. GUI] is the 
strongest, and even in the absence of strong iron blends, 
it is most often heavily blended on the blue side by 
SillljA 1892 and AlIIIA 1862 ('e.g.. lVestergaard fc Wilke^ 

fOl). The best resolved example is probably Ark 564 
eighly & Moore 2004), which indeed shows isolated 
symmetric GUI] lines at the systemic redshift, which 
mates poorly with the differential absorption hypothe- 
sis. However, this very narrow line Seyfert 1 is hardly 
representative itself. 

3) The observed excess in the red wings of the emission- 
line ratios in Fig^ is significant, and nothing similar is 
seen in Fig 13 At least for CIV, the emission peak is 
often blueshifted in BALQSOs: features in composite 
spectra of quasars with larger blueshifts correlate with 
composites of increasi ng BAL-typ c (hi towards loBAL) 
absorption (Fig. 10 in Richards et al.ll200^ This works 
against a hypothesis of differential absorption, because 
B, C, D as the more absorbed components should show 
blueshifted, not redshifted emission profiles. 

7. IMPLICATIONS AND PREDICTIONS 

An appealing aspect of our model is that it is sub- 
ject to immediate ob servational tests. X-ray observa- 
tions of BAL samples ijGreen et al.l2001trGallagher etahl 
|2002a), measure intrinsic absorbing columns of N]^*^ > 
lO^^cm"^, much higher than naively derived from the 
(typically saturated) UV BALs. Since we propose that 
A is seen via a typically less-absorbed sightline, A could 
show a less absorbed X-ray spectrum, and also larger 
fx I fopt relative to components B — D. Differential ab- 
sorption studies of lenses have been difficult in X-rays, 
due t o the close (~ 1") spacing of lensed images l)Morganl 
1200 It Ghartas 2002), but are quite feasible with wide 
lenses like SDSS J1004-I-4112. An 80ksec X-ray obser- 
vation of SDSS J1004-I-4112 was performed by Chandra 
in January 2005. The one caveat to this prediction is 
that the highly-ionized, high velocity flow proposed may 

We note that spectral differences are seen in Mgll in FiglTk. 
and Mg II absorption is the criterion for a loBAL designation. 

^ Such lines are also seen in emission in other types of AGN, but 
are most easily detec ted when narrow , for example in Narrow Line 
Seyfert Is like IZwI ILaor et aLllTOflll and quasars like 2226-3905 
JGtraham et aI.|[T99^ . 



require high S/N X-ray spectra to detect and model cor- 
rectly ijGierlihski & Dondl200l . 

Since we further propose that A is susceptible to show 
gaps in the absorber flow. X-ray measurements during 
a similar UV blue emission line asymmetry event could 
show substantially lower absorption. A very crude es- 
timate based on Figure |21 and t ypical BAL colurn ns as 
measured in the X-rays of <10^'^ (|Green et alJl200l]) sug- 
gests a column change of ~10^^. Alternatively, X-ray 
spectral fits might also reveal a change in covering factor 
of the absorption. 

Wide quasar pairs at similar redshifts are under intense 
study as possible wide lenses an d have been hunted in 
large surveys like the 2dF ( Miller _etaLll2004|) and SDSS 
flnada et al."2003l : l0gu?i et a l. 2005^. In a large enough 
statistical sample of lenses, LAMAS predicts that the 
degree of spectral difference should correlate with sepa- 
ration angle, and should be independent of proximity to 
a bright galaxy of high microlensing optical depth. The 
same can be said for the variability -f time delay interpre- 
tation, but in that case, adequate spectroscopic monitor- 
ing should reveal propagation of spectroscopic features to 
all image components. 

Perhaps the most controversial part of this proposal 
is its implication that all quasar spectra might be self- 
absorbed at some level by their smoothly outflowing 
winds, but that this fact has escaped our not ice since 
quasa rs were discovered some 40 years ago l)SchmidtJ 
1196,'^ . However, the origin of observed emission line pro- 
flies and even the quasi-power-law continuum of quasars 
remain in general poorly understood. Both their over- 
all similarity and th eir diversity are difficult to explain 
({Baldwin et alJlT995D . Absorbing winds may be a crucial 
missing component. 

Variable absorption may also be endemic. The 
salient characteristics of quasar variability seem to be 
(1) timescales of months to years (2) lower luminos- 
ity quasars have larger amplitudes of variability (3) 
at all redshifts, variability is larger at shorter wave- 
lengths (quasars are bluer when brighter) (4) variabil- 
ity amplitudes increase with redshift, indicating evolu- 
tion of the quasar po pulation or the variability mech- 
anism ()de Vries et al.l 120051. To explain (non-blazar) 
varia bility, models include accretion disk instabilities 
(e.g., iKawaguchi et al.lll99S) . so-called Poissonia n pro- 
cesse s, such as multiple supernovae (e.g., Terlevic h et alJ 
^23)1 ^iid gravitational microlensing (e.g., iHawkinsI 
199 31). Since absorption has been seen to vary sig- 
nificantly on timescales of about a month, then some 
fraction of quasar flux variability should be due to ab- 
sorbers. However, this does not require that more ab- 
sorbed QSOs are more variable. Indeed, BALQSOs show 
no different variability pr operties than QSOs generally 
l)Vanden Berk et'ani2004^ . Although the absorbers dis- 
cussed here are warm (ionized), dust may also be ex- 
pected to contribute to varying degrees which would con- 
tribute to the reddening of dimming events. According 
to Elvis c t al. ( 2002) ■ quasar outflows are natural dust 
producers. Since quasars are expected to be dustier at 
early cosmological times (jWilman et al.ll2'00ClD . this may 
naturally explain the correlation of variability with red- 
shift. 

8. CONCLUSIONS 
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We propose that small shifts in line-of-sight to quasars 
afforded by gravitational Icnsing can yield noticeable 
differences in spectroscopic line profiles, due to lat- 
eral fine structure in quasar absorbing outflows. Our 
hypothesis that lens-aided multi-angle spectroscopy of 
SDSS J1004-I-4112 probes a rotating, generally smooth, 
sheet-like flow with small-scale lateral spatial structure 
is consistent with these facts: (f ) spectroscopy of lensed 
BALQSOs generally shows similar BAL trough profiles, 
but exceptions are known (2) significant absorber column 
differences can exist between sightlines to lensed BALs 
(3) absorption columns to individual quasars are known 
to vary significantly on timescales of hours to years (4) 
one component of SDSS J1004-I-4112 shows a UV emis- 
sion line blue wing enhancement that is spectroscopically 
consistent with an unveiling event 

Models of multi-epoch spectroscopic observations of 
lensed quasar images are necessarily convolved with mod- 
els of intrinsic quasar variability, and with the as-yet 
poorly constrained structural parameters of the emit- 
ting/absorbing regions. Microlensing holds some promise 
for resolving the structure of the BELR in quasars. Un- 
fortunately, detailed and accurate lens models are re- 
quired. Worse, the size, mass, and velocity of the mi- 
crolensing caustic will remain poorly constrained because 
microlensing depends on events (caustic crossings) that 
are fundamentally unpredictable and irrcproduciblc, re- 
quiring a statistical approach to deconvolve the BELR 
structure. In contrast, if the LAMAS dominates over 
microlensing, then only the easily-measured angular sep- 
aration affects the observations. 

The most common proposed causes of spectroscopic 
differences in lens components (time variability and mi- 
crolensing) are both proven phenomena, while LAMAS 
is not. LAMAS is distinguishable in several ways from 
these, and we propose several tests above. The addition 
of LAMAS to the mix may seem to muddy the waters, 
but LAMAS should be taken seriously until disproved, 
and may in fact provide direct geometric information 
about the internal structure of AGN that has long eluded 
us. 
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Fig. 1.— SDSS J1004+4112 SPECTRA (a) At left, spectra of the 4 quasar images of SDSS J1004+4112 rescaled for 
clarity (top) reproduced from (Oguri ct al. 2004). The blue wings of emission lines are enhanced in A, most notably in CIV. 
The bottom panel highlights the differences, showing the ratio of each of B, C, and D to the A component spectrum, (b) In the 
right panel, seven epochs of the CIV emission line of SDSS J1004+4112 are reproduced directly from IRichards et al.l ll2004) . 
The spectra have been smoothed, renormalized so their peaks match, and are shown with a scaled Gaussian for reference (dotted 
line). The blue wing bump is apparent in the first 2 epochs of the A spe ctrum. Even after 21 Nov 2003, A ma i ntains a strong 
blue asymmetry. The blue bump in component A reappeared in 2004 jRichards. Johnston, fc Hennawil 120041 : IWisotzki et all 
I200I . 
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Fig. 2. — TOP: A rest-frame non-BAL composite spectrum from the SDSS from iReichard et alj[2003D (black line), along with 
composite spectra we generated by adding 25, 20, and 15% of the SDSS high-ionization BAL composite spectrum (3 spectra 
below as marked). We have taken the logarithm of flux and offset the spectra for clarity, and we highlight just the region from 
Lya to GUI], where the differences are most apparent. No BAL troughs are evident even in the 25% BAL spectrum. BOTTOM: 
The lower panel shows the ratio of these 3 part-BAL composite spectra to the non-BAL composite. While features may differ 
in detail, note the overall similarity to the spectral ratios seen in Figure^for SDSS J1004-I-4112. 



12 



Green et al. 



2366-3367™ 





n-i — 1 — 1 — 1 — I — 1 — 1— 

■ b 

■ 1 

■ L 




' ■■■!!■■■■ 1 ■■■■ 1 ' 

" 1 ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 1 1 




i" ■ 1 1 1 ■ 1 1 i ■ 1 1 1 1 i" 


SI fJ/0 W ^ c N hI 
1 ■ 


% ■ 

5 



24O0 25tW aeOO 2700 
Obaervad wovfllangth (A) 



2S00 2W>0 



Fig. 3. — 2355— 3 357 C IV region of the spectrum of the quasar PKS J2355-3357, with multi-component fits from 
IKuraszkiewicz et alJ (12003) . The blue wings of CIV emission line is enhanced here, with a profile similar to that in 
SDSSJ1004+4112. 
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Fig. 4. absor ber MODEL. LEFT: A cartoon of a stratified wind outflow model JMurrav fc Chiansll997llElvisl2000l: 
lEverett et al.ll2nn2ri shows 2 observers at slightly different line-of-sight angles to the outflow. In these models, wind-embedded 
clouds/streams act as emitters (abso rbers) when we look at (through) them. RIGHT: A Gaussian emission line absorbed by 
a standard BAL outflow model (e.g.. IAravlll999D along 2 slightly different sightlines as shown at left. 



